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Abstract. 

We study the pressure and temperature dependences of the dielectric relaxation 
of two molecular glassforming liquids, dibutyl phtalate and m-toluidine. We focus 
on two characteristics of the slowing down of relaxation, the fragility associated with 
the temperature dependence and the stretching characterizing the relaxation function. 
We combine our data with data from the literature to revisit the proposed correlation 
between these two quantities. We do this in light of constraints that we suggest to 
put on the search for empirical correlations among properties of glassf ormers. In 
particular, argue that a meaningful correlation is to be looked for between stretching 
and isochoric fragility, as both seem to be constant under isochronic conditions and 
thereby reflect the intrinsic effect of temperature. 
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1. Introduction 

With the goal of better understanding the physics of glasses and of glass formation, 
there has been a continuing search for empirical correlations among various aspects 
of the phenomenology of glassformers. The most distinctive feature of glass 
formation being the rapid increase with decreasing temperature of the viscosity 
and relaxation times, correlations have essentially been looked for between the 
characteristics of the latter and other thermodynamic or dynamic quantities. Angell 
coined the term "fragility" to describe the non-Arrhenius temperature dependence 
of the viscosity or (alpha) relaxation time and the associated change of slope on an 
Arrhenius plot He noticed the correlation between fragility and amplitude of 
the heat-capacity jump at the glass transition. Earlier, the Adam-Gibbs approach 
was a way to rationalize the correlation between the viscosity increase and the 
configurational or excess entropy decrease as one lowers the temperature 0. 
Since then, a large number of empirical correlations between "fragility" and other 
properties of the liquid or of the glass have been found: for instance, larger fragility 
{i.e., stronger deviation from Arrhenius behavior) has been associated with (i) a 
stronger deviation of the relaxation functions from an exponential dependence on 
time (a more important "stretching") IfSl, (ii) a lower relative intensity of the boson 
peak |4l, (iii) a larger mean square displacement at Tg |[5l, (iv) a smaller ratio of 
elastic to inelastic signal in the X-ray Brillouin-spectra lHJ, (v) a larger Poisson ratio 
|[7| and (vi) a stronger temperature dependence of the elastic shear modulus. Goo, in 
the viscous liquid ||8]]. 

For useful as they may be to put constraints on proposed models and theories of 
the glass transition, such correlations can also be misleading by suggesting causality 
relations where there are no such things. It seems therefore important to assess the 
robustness of empirically established correlations. In this respect, we would like to 
emphasize a number of points that are most often overlooked: 

1) Fragility involves a variation with temperature that a priori depends on the 
thermodynamic path chosen, namely constant pressure (isobaric) versus constant 
density (isochoric) conditions. On the other hand, many quantities that have been 
correlated to fragility only depend on the thermodynamic state at which they are 
considered: this is not the case for the variation of the excess entropy or of the 
shear modulus, nor for the jump in heat capacity measured in differential scanning 
calorimetry, which are all path dependent, but the other properties are measured 
either at Tg, the glass-transition temperature, or in the glass, where they also relate 
to properties of the liquid as it falls out of equilibrium at Tg (there may be a residual 
path dependence due to the nonequilibrium nature of the glass, but it is quite 
different from that occuring in the liquid). Which fragility then, isobaric or isochoric, 
should best be used in searching for correlations ? 

2) The quantities entering in the proposed correlations are virtually always 
considered at Tg. This is the case for the commonly used measure of fragility, the 
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"steepness index", which is defined as the slope of the temperature dependence of 
the alpha-relaxation time on an Arrhenius plot with T scaled by Tg |9l- Tg is of 
course only operationally defined as the point at which the alpha-relaxation time (or 
the viscosity) reaches a given value, say 100 seconds for dielectric relaxation. The 
correlated properties are thus considered at a given relaxation time or viscosity. What 
is the fate of the proposed correlations when one studies a different value of the relaxation time 

3) Almost invariably, comparisons involve properties measured at atmospheric 
pressure, for which the largest amount of data is available. Since, as discussed in 
the preceding point, the properties are also considered at a given relaxation time, an 
obvious generalization consists in studying the validity of the reported correlations 
under "isochronic" {i.e., constant relaxation time ) conditions, by varying the control 
parameters such that the relaxation time stays constant. How robust are then the 
correlations when one varies, say, the pressure along an isochrone ? In light of the above, 
our contention is that any putative correlation between fragility and another property 
should be tested, as far as possible, by varying the reference relaxation time, by 
varying the thermodynamic state along a given isochrone, and by changing the 
thermodynamic path along which variations, such as that defining the fragility, are 
measured 

A better solution would certainly be to correlate "intrinsic" properties of 
glassformers that do not depend on the chosen state point or relaxation time. A step 
toward defining such an "intrinsic" fragility has been made when it was realized 
that the temperature and the density dependences of the alpha-relaxation time and 
viscosity of a given liquid could be reduced to the dependence on a single scaling 
variable, X = e{p)/T, with e{p) an effective activation energy characteristic of the 
high- temperature liquid IITOl [TT| . Evidence is merely empirical and is supported 
by the work of several groups for a variety of glassforming liquids and polymers 
(TOl [III [121 [131 [11 [151 [111. The direct consequence of this finding is that the fragility 
of a liquid defined along an isochoric path is independent of density: the isochoric 
fragility is thus an intrinsic property, contrary to the isobaric fragility. Although one 
could devise ways to characterize the isochoric fragility in a truly intrinsic manner, 
independently of the relaxation time, the common measure through the steepness 
index (see above) still depends on the chosen isochrone. In looking for meaningful 
correlations to this isochoric steepness index, it is clear however that one should 
discard quantities that vary with pressure (or equivalently with temperature) under 
isochronic conditions. As we further elaborate in this article, the stretching parameter 
characterizing the shape of the relaxation function (or spectrum) is a priori a valid 
candidate, as there is some experimental evidence that it does not vary with pressure 
along isochrones | [T7| . 

The aim of the present work is to use the knowledge about pressure and 
temperature dependences of the liquid dynamics to test the robustness of proposed 
correlations between fragility and other properties. This is a continuation of the work 
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presented in reference IITSB , where the focus was mainly on correlations between 
fragility of the liquid and properties of the associated glass. In this paper we 
specifically consider the correlation between fragility and stretching. The reported 
correlation between the two is indeed one of the bases of the common belief that both 
fragility and stretching are signatures of the cooperativity of the liquid dynamics. 

We present new dielectric spectroscopy data on the pressure dependence of 
the alpha relaxation of two molecular glassforming liquids, dibutyl phthalate (DBF) 
and m-toluidine. We express the alpha-relaxation time as a function of the scaling 
variable X = e{p)/T and evaluate the density dependence of e{p) as well as the 
isochoric fragility. We also study the spectral shape and its pressure dependence 
along isochronic lines. We spend some time discussing the methodological aspects 
of the evaluation of the fragility and of the stretching from experimental data, as well 
as that of the conversion from P, T to P, p data. This provides an estimate of the error 
bars that one should consider when studying correlations. Finally, by combining 
our data with literature data we discuss the robustness of the correlation between 
fragility and stretching along the lines sketched above. 

The paper is structured as follows. Section |2l introduces some concepts and 
earlier developments that are central for the discussion. In section |3] we present 
the experimental technique. Section |4] is devoted to the pressure, temperature and 
density dependence of the relaxation time. In section|5]we analyze the spectral shape 
and its pressure and temperature dependence. Finally, in section [6] we combine 
the current results with literature data to assess the relation between fragility and 
stretching, stressing the need to disentangle temperature and density effects. Two 
appendices discuss some methodological points. 

2. Background 

2.1. Isochoric and isobaric fragilities 

The fragility is a measure of how much the temperature dependence of the alpha- 
relaxation time (or alternatively the shear viscosity) deviates from an Arrhenius form 
as the liquid approaches the glass transition. The most commonly used criterion is 
the so called steepness index. 



where the derivative is evaluated at Tg and is expressed in seconds. 
Conventionally, the liquid is referred to as strong if m is small, that is 17 — 30, 
and fragile if m is large, meaning roughly above 60. In the original classification 
of fragility it was implicitly assumed that the relaxation time (or viscosity) was 
monitored at constant (atmospheric) pressure, as this is how the vast majority of 
experiments are performed. The conventional fragility is therefore the (atmospheric 
pressure) isobaric fragility, and, as indicated in Eq. [U, the associated steepness index 




(1) 
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is evaluated at constant pressure. However, the relaxation time can also be measured 
as a function of temperature along other isobars, and this will generally lead to 
a change in ntp. Moreover, it is possible to define an isochoric fragility and the 
associated index, ntp, obtained by taking the derivative at constant volume rather 
than at constant pressure. The two fragilities are straightforwardly related via the 
chain rule of differentiation, 

dp 



51ogio(T«) 
mp = rup + - " 



dp 



{T 

P 



dTg/T 

when both are evaluated at the same point {Tg{P), p{P,Tg{P))). The isochoric 
fragility, ntp, describes the intrinsic effect of temperature, while the second term on 
the right hand side incorporates the effect due to the change of density driven by 
temperature under isobaric conditions. It can be shown that the above relation can 
be rewritten as 

nip = nip(l — ap/ocr) (2) 

where the unconventional is the isochronic expansivity IIT9B , i.e., the expansivity 
along a line of constant alpha-relaxation time Ta (the Tg line being a specific 
isochrone). The above result is purely formal and contains no assumptions. The 
implication of the result is that ntp is larger than ntp ii ap > and Ur < 0. It is well 
known that > in general. The fact that cx-r is negative arises from the empirical 
result that the liquid volume always decreases when heating while following an 
isochrone. 

Within the last decade a substantial amount of relaxation-time and viscosity 
data has been collected at different temperatures and pressures /densities. On the 
basis of the existing data, it is reasonably well established that the temperature and 
density dependences of the alpha-relaxation time can be expressed in a scaling form 

as cniiniiiaiisiiiiiisEi- 

r«(p,r)=f(^). (3) 

It is seen directly from Eq. |3]that X(p, T) = e{p)/T, when evaluated at Tg, has the 
same value at all densities (X^^ = e{p)/Tg{p)) if Tg{p) is defined as the temperature 
where the relaxation time has a given value (e.g., = 100 s). Exploiting this fact, 
it is easy to show [TTl |20ll that the scaling law implies that the isochoric fragility is 
independent of density. For instance, the isochoric steepness index, when evaluated 
at a Tg corresponding to a fixed relaxation time, is given by 



^ _ dlogio(T«; 



'9 



dTg/T 



[T = Tg) = F'iXg)^{T = Tg) = XgF'iXg).{A) 



p 

The fact that the relaxation time is constant when X is constant means that the 
isochronic expansion coefficient is equal to the expansion coefficient at constant 
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X. Using this and the general result (^|^) ^ (^||) ^ (^g) = -1, it follows that 



= _j dloge(p) 
ar ^ d log p 

which inserted in Eq. |2] leads to 



Tg , \ (5) 



, dlog£'(p)' 
nip = nip [1 + apTg—f , (6) 



dlog p 

where mp, rup and ap are evaluated at Tg. 

When liquids have different isobaric fragilities, it can be thought of as due to 
two reasons: a difference in the intrinsic isochoric fragility, ntp, or a difference in the 
relative effect of density, characterized by apTg and the parameter x = ^d'log ^ • 
analyze the data in this frame. 

2.2. Relaxation-time dependent fragility 

The following considerations hold for isochoric and isobaric conditions alike. The p 
or P subscript are therefore omitted in this section. 

The fragility is usually characterized by a criterion evaluated at Tg, i. e., the 
temperature at which the relaxation time reaches Ta;=100 s-1000 s. The same criterion, 
e.g. the steepness index, can however equally well be evaluated at a temperature 
corresponding to another relaxation time, and this is also found more often in 
literature, mainly to avoid the extrapolation to long times. So defined, the "fragility" 
for a given system can be considered as a quantity which is dependent of the 
relaxation time at which it is evaluated: 

m(r) = ^!§^(T = T.) (7) 

where Tq;(Tt-) = t defines the temperature T^. {Tg is a special case with t ^100 
s-lOOOs .) 

An (extreme) strong system is a system for which the relaxation time has an 
Arrhenius behavior, 

E 



T«(r) = Too exp y-jrj ' (8) 

where Eoo is a temperature and density independent activation energy (measured 
in units of temperature). Inserting this in the expression for the relaxation-time 
dependent steepness index (Eq. [Z]) gives 

nistrong{r) = logio (t/Too) (9) 

which gives the value mstwngi"^ = 100s) = 15 (assuming Zogio (Too/sec) = —13) and 
decreases to nistrongi 

T — Too) — as the relaxation time is decreased. This means that 
even for a strong system the steepness index is relaxation-time dependent. In order 
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to get a proper measure of departure from Arrhenius behavior it could therefore be 
more adequate to use the steepness index normalized by that of a strong system: 



mJr) = ^(^^ = ^^-/^ (10) 

nistrongir) logio (t/Too) 

nin{r) will take the value 1 at all relaxation times in a system where the relaxation 
time has an Arrhenius behavior. Such a normalized measure of fragility has been 
suggested before | |2ni22ll23| . For instance, Olsen and coworkers ||23B have introduced 
the index 



where E(T) is a temperature dependent activation energy defined by E(r) = 
Tin (tcc/ Tin fty). The relation between the steepness index and the Olsen index is 

m 

J(t) = ^^J^ - 1 = m„(T) - 1 (12) 

I{r) takes the value for strong systems at all relaxation times. Typical glass 
forming liquids display an approximate Arrhenius behavior at high temperatures 
and short relaxation times; in this limit I(t) = and it increases as the temperature 
dependence starts departing from the Arrhenius behavior. Typical values of I at 
Tg{r = 100s) are ranging from 1=3 to 1=8, corresponding to steepness indices of 
m=47 to m=127. 

Finally, we note in passing that relaxation-time independent measures of 
fragility can be formulated through fitting formulae: this is the case for instance of 
the fragility parameter D in the Vogel-Tammann-Fulcher (VTF) formula or of the 
frustration parameter B in the frustration-limited domain theory [[24||. 



3. Experimentals 

The dielectric cell is composed of two gold-coated electrodes separated by small 
Teflon spacers. The distance between the spacers is 0.3 mm and the area is 5.44 
cm^ giving an empty capacitance of 16 pF. The electrodes are totally immersed in 
the liquid sample, which is sealed from the outside by a Teflon cell. The electric 
contacts are pinched through the Teflon. The compression is performed using liquid 
pentane, which surrounds the Teflon cell from all sides. The Teflon cell has one end 
with a thickness of 0.5 mm in order to insure that the pressure is well transmitted 
from the pentane to the liquid sample. The pressure is measured by using a strain 
gauge. The cooling is performed by a flow of thermostated cooling liquid running 



On the correlation between fragility and stretching in glassforming liquids 



8 



inside the autoclave. The temperature and the temperature stability are monitored 
by two PTIOO sensors placed 2cm and 0.3 cm from the sample. The temperature is 
held stable within ±0.1 degree for a given isotherm. The temperature during the 
time it takes to record a spectrum is stable within ±0.01 degree. 

The setup insures a hydrostatic pressure because the sample is compressed 
from all sides. It is moreover possible to take spectra both under compression and 
decompression. By doing so and returning to the same P — T condition after several 
different thermodynamic paths, we have verified that there was no hysteresis in the 
pressure dependence of the dynamics. This serves to confirm that the liquid is kept 
at thermodynamic equilibrium at all stages. 

The capacitance was measured using a HP 4284A LCR-meter which covers the 
frequency range from 100 Hz to 1 MHz. The low-frequency range from 100 Hz to 1 
Hz is covered using a SR830 lockin. 

The samples, dibutyl phthalate (DBF) and m-toluidine, were acquired from 
Sigma- Aldrich. The m-toluidine was twice distilled before usage. The DBF was used 
as acquired. 

Liquid m-toluidine was measured on one isotherm at 216.4 K. DBF was 
measured along 4 different isotherms, 205.5 K, 219.3 K, 236.3 K and 253.9 K, at 
pressures up to 4 kbar. DBF was moreover measured at different temperatures along 
two isobars: atmospheric pressure and 230 MFa. The pressure was continuously 
adjusted in order to compensate for the decrease of pressure which follows from 
the contraction of the sample due to decreasing temperature. It is of course always 
possible to reconstruct isobars based on experiments performed under isotherm 
conditions. However, such a procedure mostly involves interpolation of the data, 
which is avoided by performing a strictly isobaric measurement. For DBF we have 
obtained relaxation-time data at times shorter than 10~^'^ s by using the high- 
frequency part of the spectrum and assuming time-temperature and time-pressure 
superposition (TTFS). Although TTFS is not followed to a high precision (see section 
15. 1|) , the discrepancies lead to no significant error on the determination of the 
relaxation time. This is verified by comparison to atmospheric-pressure data from 
the literature (see figured]). 

4. Alpha-relaxation time and fragility 

4.1. Dibutyl phtalate 

The DBF data at atmospheric pressure is shown in figure [1] along with literature 
results. Tg{Patm) = 177 K, when defined as the temperature at which = 100 
s. We also present the data taken at P = 230 MFa in this figure. It is clearly seen that 
Tg increases with pressure. An extrapolation of the data to t^, = 100 s gives Tg = 200 
K for P = 230 MFa, corresponding to dTg/dP ~ 0.1 KMFa~^. This corresponds well 
to the pressure dependence of Tg (at = 1 s) reported by Sekula et al. |[25l , based on 
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measurements taken at pressures higher than 600 MPa. The dependence is however 
stronger than that reported by Fujimori et al. 1261 based on isothermal calorimetry, for 
which dTg/dP = 0.06 KMPa~^. This indicates that the calorimetric and the dielectric 
relaxations may have somewhat different dependences on pressure. 

In figure l2l we illustrate the determination of Tg and of the steepness index mp 
for the atmospheric-pressure data, using the part of the data of figure [1] whith a 
relaxation time longer than a millisecond. Along with the data we show the VTF 
fit from Sekula et al. [|25| extrapolated to low temperatures, which gives Tg = 177 A K 
and mp = 84. We have also performed a new VTF fit restricted to the data in the 10~^ 
s— 10^ s region. The result of this fit yields Tg = 176.1 K and mp = 79. Finally, we 
have made a simple linear estimate of logiQXoc a function of l/T in the temperature 
range shown in the figure. This linear slope fits the data close to Tg better than any of 
the VTF fits. The corresponding glass transition temperature and steepness index are 
Tg = 176 K and mp = 65. This illustrates that the determination of Tg is rather robust 
while this is less so for the steepness index. This latter depends on how it is obtained, 
and the use of extrapolated VTF fits can lead to an overestimation. (Of course, a VTF 
fit made over a very narrow range, e.g. lO'h - 10^ s, will agree with the linear fit, 
because it becomes essentially linear over the restricted range.) The fragility of DBF 
has earlier been reported to be mp = 69 BH, based on the data of Dixon et al. [|27|. We 
take mp = 67 as a representative value. 

The relaxation-time data along four different isotherms are displayed as a 
function of pressure in figure IH 

In order to separate the relative effects of density and temperature it is 
convenient to express the relaxation time as a function of density and temperature 
rather than pressure and temperature. To do this, we need the pressure and 
temperature dependences of the density. However, for liquid DBF such data is only 
available at high temperature ||28| . 

In order to extrapolate the equation of state to low temperature we have applied 
the following scheme. When calculated from the data in Ref. ||28|, the expansion 
coefficient ap shows a weak decrease with decreasing temperature. We therefore 
assume that the temperature dependence of ap is linear over the whole temperature 
range and integrate with respect to temperature to obtain the density along the 
atmospheric-pressure isobar. In the whole temperature range of Ref. ||28l , the 
pressure dependence of the density is well described by fits to the Tait equation 
with temperature-dependent adjustable parameters "c" and "h" l|29ll (which are 
directly related to the compressibility and its first-order pressure derivative). We 
have linearly extrapolated the temperature dependence of these parameters and 
used the Tait equation to calculate the pressure dependence along each isotherm. 
Extrapolating the derivatives rather than the density itself is expected to lead to 
smaller errors on the latter. In addition, we have checked that this procedure gives 
physically reasonable pressure and temperature dependences of the expansivity and 
of the compressibility Il30ll . 
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Figure H] shows the density dependence of the alpha-relaxation time along the 
four different isotherms, the atmospheric-pressure isobar and the 230 MPa isobar. 
We have also included the room-temperature dielectric data of Paluch et al. ||3T|. 
For DBF the viscosity data and the dielectric relaxation time do not decouple under 
pressure ||25|, and we have therefore also included the room-temperature viscosity 
data of Cook et al. ||3T1 . 

In figure |5] we show the data of figure H] plotted as a function of the scaling 
variable p^/T, choosing for x the value that gives the best collapse for the data of this 
work. This corresponds to testing the scaling in equation |3]by assuming that e(p) is a 
power law. The data taken at low density collapse quite well with x = 2.5, while this 
is not true for the data of Faluch 1311 taken at densities higher than approximately 1.2 
g/ cm^. It is possible to make all the data collapse by allowing e{p) to have a stronger 
density dependence at higher densities. In figure [6] we show the data as a function 
of e(p)/T, where we have constructed the density dependence of e(p) in order to get 
a good overlap of all the data (we did not look for the best collapse, but merely 
evaluated the change of the isochronic expansivity: see section The resulting 
density dependence of e(p) is shown in figure [6] along with the p^-^ power law. Note 
that the quality of the data collapse depends only on the density dependence of 
e(p) not on its absolute value. The constructed e{p) has an apparent "power-law" 
exponent x(p) = dloge{p)/dlogp that increases from 1.5 to 3.5 with density in the 
range considered. In any case, the absence of collapse in figure |5]cannot be explained 
by errors in estimating the FVT data: this is discussed in more detail in Appendix 



Appendix A 



As a last note regarding the e(p)/T-scaling in figure [6l we want to stress that we 
cannot test the scaling (Eq. [S]) in the density range above 1.25g/cm'^ where there is 
only one set of data. (This is why we did not attempt to fine tune e{p) to find the best 
collapse, see above.) Indeed, with a unique set of data in a given range of density it 
is always possible to construct e{p) in this range to make the data overlap with data 
taken in other density ranges. 

We have determined the ratio between the isochoric fragility and the isobaric 
fragility at atmospheric pressure by calculating ar along the isochrone of 100 s and 
inserting it in Eq. |2l This leads to mp/ntp ~ 1.2, when ntp. In figure [71 we show the 
isobaric data taken at atmospheric pressure and at 230MFa scaled by their respective 
Tg{P). No significant pressure dependence of the isobaric fragility is observed when 
going from atmospheric pressure to 230 MFa, which is consistent with the result of 
reference ||25|. The pressure independence of mp is connected to the relatively low 
value of mp/rup = 1.2 (typical values are 1.1-2 ||TT|): mp is pressure independent and 
the ratio ntp / nip cannot be lower than one (see Eq. |2]), so that mp can at most decrease 
by 20% from its atmospheric-pressure value. Such a change would almost be within 
the errorbar of the determination of mp from the data at 230MFa (see the discussion 
earlier in this section). 
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Figure 1. Temperature dependence of the alpha-relaxation time (from dielectric 
measurements, Ta = 1 /cUpeak) of liquid DBF at atmospheric pressure and at 230 MPa 
(Arrhenius plot). Data at atmospheric pressure from other groups are also included: 
unpublished data from Nielsen Il32l , the VTF fit of Il25ll shown in the range where it 
can be considered as an interpolation of the original data, and data taken from figure 
2(a) in reference [27]. 
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Figure 2. Atmospheric-pressure data of figure [T] with relaxation times longer than a 
millisecond (symbols). Also shown are the VTF fit from reference ||25l extrapolated to 
low temperatures (dashed-dotted line), a new VTF fit made by using data in the 10~^ 
s — 10^ region (dashed line), and estimated slope of the data in the long-time region 
(full line). The Tg's estimated from these three methods are very similar, whereas the 
fragility varies significantly from m = 65 to m = 85. 
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Figure 3. Alpha-relaxation time of DBP (from dielectric measurements, Xa = 1 /cupeak) 
as a function of pressure along 4 different isotherms (log-linear plot). 
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Figure 4. Logarithm of the alpha-relaxation time of DBP versus density (see the text 
regarding the calculation of density). Included are data from this work along with 
dielectric data from figure 3 in reference IISTl , and viscosity data from reference |29l . 
The viscosity data are shifted arbitrarily on the logarithmic scale in order to make the 
absolute values correspond to the dielectric data of reference ISH , which are taken at 
the same temperature. 
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Figure 5. The alpha-relaxation times shown in figurelUplotted as a function of p^'^/T. 
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b) 




Figure 6. (a) The alpha-relaxation times shown in figure |4] plotted as a function of 
X = e{p)/T, with increasing dloge{p) /dlogp as p increases, (b) Density-dependent 
activation energy e(p) (dashed line) used in the scaling variable X = e{p)/T for 
collapsing data in (a) (the associated x(p) = dloge{p) /dlogp increases from 1.5 to 
3.5 in the density range under study). We also display the power law giving the best 
scaling, p^-^, at low density (full line). 
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Figure 7. Arrhenius plot of the alpha-relaxation time of DBP at atmospheric pressure 
and at 230 MPa, when the temperature is scaled with the pressure dependent Tg, 
Tg{Patm) = 176K and Tg(230MPa)=200K. As in figure [H data from other groups are 
also included: unpublished data from Nielsen Il32l , the VTF fit of l25ll shown in the 
range where it can be considered as an interpolation of the original data, and data 
taken from figure 2 (a) in reference ||27). 



4.2. m-Toluidine 

The glass transition temperature at atmospheric pressure is = 187 K (for t^; = 100 
s) and the isobaric fragility based on dielectric spectra is reported to be mp = 82 ± 3 
||33l |34B . (There has been some controversy about the dielectric relaxation in m- 



toluidine, see reference p3|l and references therein.) 

In the inset of figure |8] we show the pressure-dependent alpha-relaxation time at 
216.4 K. Extrapolating the data to t« = 100 s leads to = 340 ± 10 MPa, which is 
in agreement with the slope, dTg/dP = 0.45 KMaP~^, reported for the calorimetric 
glass transition in ||35| . This indicates that the decoupling between the timescales of 
dipole relaxation and of calorimetric relaxation which appears under pressure in the 
case of DBP is not present in m-toluidine in this pressure range. 

As for DBP, we wish to convert the temperature and pressure dependences of 
the relaxation time to the temperature and density dependences. Density data are 
available along four isotherms in the 278.4 K —305.4 K range for pressures up to 
300 Mpa ||36|. Tait fits and thermal expansivity in this range were extrapolated by 
using the scheme described above for DBP in order to determine density both as a 
function of temperature down to Tg, and as a function of pressure on the 216.4 K 
isotherm. In figure |8] we show the alpha-relaxation time as a function of density. The 
data taken at atmospheric pressure and the data taken along the 216.4 K isotherm 
cover two different ranges in density. It is therefore not possible from this data to 
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verify the validity of the scaling in X = e{p)/T. We therefore assume that the scaling 
is possible. Moreover, due to the paucity of the data we describe e{p) by a simple 
power law, e{p) = p^. We find the exponent x by exploiting the fact that the scaling 
variable X = e(p)/r is uniquely fixed by the value of the relaxation time; applying 
this at Tg, namely setting Xg{Patm) =X^(216K), leads to x = 2.3 and gives a ratio of 
mp/nip = 1.2. 




Figure 8. Logarithm of the alpha-relaxation time of m-toluidine as a fimction of 
density along the isotherm T = 216.4K (symbols). The VTF fit of the atmospheric- 
pressure data of reference 1331 is also shown in the range where the fit can be 
considered as an interpolation of the data (dashed line). The inset shows the alpha- 
relaxation time of m-toluidine as a function of pressure along the isotherm T=216.4 
K. 



5. Spectral shape and stretching 

The shape of the relaxation function (or spectrum), most specifically its distinctly 
nonexponential (or non-Debye) character in the viscous regime, is taken as one of the 
important features of glassforming materials. Characterizing and quantifying this 
effect is however not fully straightforward and has led to diverging interpretations. 
First of all, the shape of the relaxation function or spectrum may change with the 
experimental probe considered. Even when restricting comparison to a single probe, 
here, dielectric relaxation, there is no consensus on how to best characterize the 



shape. We discuss in appendix Appendix B various procedures that are commonly 



used and we test their validity on one representative spectrum. For reasons detailed 
in that appendix, we focus in the following on the Cole-Davidson fitting form. 
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5.1. Dibutyl phtalate 

The frequency-dependent dielectric loss for a selected set of different pressures 
and temperatures is shown in figure HI The first observation is that cooling and 
compressing have a similar effect as both slow down the alpha relaxation and 
separate the alpha relaxation from higher-frequency beta processes. The data 
displayed are chosen so that different combinations of temperature and pressure 
give almost the same relaxation time. However, the correspondence is not perfect. 
In figure [lO] we have thus slightly shifted the data, by at most 0.2 decade, in order 
to make the peak positions overlap precisely. This allows us to compare the spectral 
shapes directly. It can be seen from the figure that the shape of the alpha peak itself is 
independent of pressure and temperature for a given value of the alpha-relaxation 
time {i.e., of the frequency of the peak maximum), while this is not true for the 
high-frequency part of the spectra which is strongly influenced by the beta-relaxation 
peak (or high-frequency wing). When comparing datasets that have the same alpha- 
relaxation time one finds that the high-frequency intensity is higher for the pressure- 
temperature combination corresponding to high pressure and high temperature. 

In figure [H] we show all the datasets of figure |9] superimposed and we zoom 
in on the region of the peak maximum. The overall shape of the alpha relaxation is 
very similar at all pressures and temperatures. However, looking at the data in more 
detail, one finds a significantly larger degree of collapse between spectra which have 
the same relaxation time, whereas a small broadening of the alpha peak is visible as 
the relaxation time is increased. At long relaxation times there is a perfect overlap of 
the alpha-relaxation peaks which have the same relaxation time. At shorter relaxation 
time, logio{uJmax) ~ 5, the collapse is not as good: the peak gets slightly broader 
when pressure and temperature are increased along the isochrone. In all cases, the 
alpha peak is well described by a Cole-Davidson (CD) shape. The ISqd goes from 
0.49 to 0.45 on the isochrone with shortest relaxation time and decreases to about 0.44 
close to Tg at all pressures. On the other hand, a Kolraush-William-Watts (KWW) fit 
close to Tg gives I3kww = 0.65. A detailed discussion of the fitting procedures and of 



the relation between CD and KWWW descriptions is given in appendix Appendix B 



5.2. m-toluidine 

The frequency-dependent dielectric loss of m-toluidine for several pressures along 
the T=216.4 K isotherm is shown in figure [121 The data are then superimposed by 
scaling the intensity and the frequency by the intensity and the frequency of the peak 
maximum, respectively: this is displayed in figure [13l When zooming in (figure [13] 
(b) we still see almost no variation of the peak shape. For the present set of data, 
pressure-time-superposition is thus obeyed to a higher degree than in DBF, and the 
changes are too small to give any pressure dependence in the parameters when fitting 
the spectra. The Cole-Davidson fit to the m-toluidine gives ^CD = 0-42 (see also 



appendix Appendix Bl. Mandanici |[33[1 and coworkers have reported a temperature 
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Figure 9. Log-log plot of the frequency-dependent dielectric loss of DBF. The curves 
can be sorted in 4 groups, each group having roughly the same peak frequency; from 
right to left: (i) Red dashed-dotted curve: T=253.9 K P=320 MPa; black dots: T= 236.3 
K and, from right to left, P=153 MPa, P=251 MPa, P=389 MPa; full blue line: T=219.3 
K and, from right to left, P=0 MPa, P=108 MPa, P=200 MPa, P=392 MPa; magenta 
dashed curve: T=206 K and, from right to left, P=0 MPa, P=85 MPa, P=206 MPa. 
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Figure 10. Same dielectric loss data of DBP as in figure |9] with a slight shift of the 
peak frequencies (less than 0.2 decade) to make the data taken under quasi isochronic 
conditions precisely coincide. The symbols are the same as in figurelH but the data at 
T=206 K and P=206 MPa and 219.3 K and P=392 MPa are not shown. 



independent value of /3cd = 0.45 for data taken at atmospheric pressure in the 
temperature range 190 K-215 K, a value that is compatible with ours. 

6. Discussion 

6.1. Correlations with fragility 

As discussed in the Introduction, the temperature dependence of the alpha-relaxation 
time (or of the viscosity) is usually considered as the most important phenomenon 
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Figure 11. Same dielectric -loss data as in figures |9] and [TOl with the frequency and 
intensity now scaled by the values at the maximum. We show only a 1.5 decade 
in frequency in order to magnify the details. Notice a small broadening as the 
characteristic relaxation time increases: Blue dashed-dotted line are three different 
data sets with logio^max « 2.6 (P=320 MPa,T=253.9 K and P=153 MPa,T=236.3 K 
and P=0 Mpa,T=219.3 K). Red full lines are three data sets with logiQ-Vmax ~ 4.1 
(P=251 MPa,T=236.3 K and P=108 MPa,T=219.3 K and P=0 Mpa,T=205.6 K). Green 
dashed lines are three data sets with logiQ-Vmax ~ 5.2 (P=339 MPa,T=236.3 K and 
P=200 MPa,T=219.3 K and P=85 Mpa,T=205.6 K). 




Figure 12. Log-log plot of the frequency-dependent dielectric -loss of m-toluidine at 
T=216.4K and pressures MPa, 59 MPa, 79 MPa, 105 MPa, 122 MPa, 142 MPa, 173 
MPa and 191 MPa. The peak shifts left as pressure is applied. Lines are guides to the 
eye. 



to understand in glass science. Isobaric fragility is then often used to characterize 
the viscous slowing down and its measures, such as the steepness index, are 
then considered as fundamental parameters. Many studies have been aimed at 
investigating which other properties of the liquid and of the associated glass correlate 
to fragility. Such correlations have been empirically established by comparing rather 
large sets of systems covering a wide spectrum of fragilities. 

In the literature, the finding of a correlation between fragility and some other 
property is always interpreted as indicating that the property in question is related 
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Figure 13. Same dielectric-loss data as in figure [121 now with the intensity and the 
frequency scaled by the values of the peak maximum. Figure (b) shows a zooming in 
of the data in (a) to focus on the alpha-relaxation region near the peak maximum. 



to the effect of temperature on the structural relaxation. However, when cooling a 
liquid isobarically two effects contribute to the slowing down of the dynamics: the 
decrease of temperature and the associated increase of density. Hence, the isobaric 
fragility is a combined measure of the two effects. It is of course the underlying goal 
that the proposed correlations be used as guidelines and tests in the development of 
theories and models for the glass transition. It is therefore important to clarify if the 
correlations result from, and consequently unveil information on, the intrinsic effect 
of temperature on the relaxation time, the effect of density, or a balanced combination 
of the two. 

Eq.s [6] and |2] show how isobaric fragility can be decomposed into two 
contributions, that of temperature being given by nip and the relative effect of density 
on relaxation time characterized by <^pTg^^^^^. Isobaric measurements do not give 

access to ntp nor to <^pTg^^^^^^ independently, but the relevant information can be 
obtained from data taken under pressure, as we have shown for the data presented 
here. From this information it becomes possible to revisit the correlations between 
fragility and other properties |T8|. The underlying idea is that a property supposed to 
correlate to the effect of temperature on the relaxation time should more specifically 
correlate to the isochoric fragility, rup, than to the isobaric one, rtip. 

As also stressed in the Introduction, it is instructive to consider the evolution 
of the empirically established correlations with pressure. As shown in section 12.11 
rup is constant, i.e., is independent of density and pressure, when it is evaluated 
at a pressure (or density) dependent Tg corresponding to a given relaxation time. 
Nonetheless, it follows from Eq. [6l that the isobaric fragility will in general change 
due to the pressure dependence of ocpTg '^^diogp^ ■ increases with pressure, ocpTg{P) 

decreases, whereas '^^Ifg^p^ = x is often to a good approximation constant (the 
DBF case at high pressure discussed in section 14.11 is one exception). As a result, 
the pressure dependence of mp is nontrivial. DBF, which we have studied here, 
shows no significant pressure dependence of the isobaric fragility, while the general 
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behavior seen from the data compiled by Roland et al. ||T3l is that the isobaric fragility 
decreases or stays constant with pressure, with few exceptions. This seems to indicate 
that the decrease of apTg{P) usually dominates over the other factors. 

The properties that are correlated to fragility will a priori also depend on pressure 
or density. However if a property is related to the pure effect of temperature on 
the relaxation time, and therefore correlates to tUp, then it should be independent of 
density when evaluated along an isochrone (usually the glass transition line Tg), as 
tUp itself does not depend on density. 



6.2. Stretching and fragility 

One of the properties that has been suggested to correlate to the fragility is the 
nonexponential character of the relaxation function, usually expressed in terms of 
the stretching parameter I3kww- 

The data we have reported here confirm the earlier finding | [T7| that the spectral 
shape of the alpha relaxation does not vary when pressure is increased while keeping 
the relaxation time constant. This leads us to suggest that, if a correlation between 
fragility and stretching does exist, this latter should better correlate to the isochoric 
fragility which is also independent of pressure than to the isobaric fragility. To test 
this hypothesis we have collected data from literature reporting isobaric fragility and 
stretching of the relaxation at Tg. We consider here a description of the shape of the 
relaxation function in terms of the KWW stretching parameter ISfcww- This choice 
is made because it is convenient to use a characterization with only one parameter 



for the shape (see appendix Appendix B for a discussion and the connection with the 
Cole-Davidson description used above) and because iS^ww is the most reported of 
the liquids where ntp is also available. The compilation of this data is shown in table 
I and in figures [T4l and [T5] where both the isobaric fragility at atmospheric pressure 
(Fig. and isochoric fragility (Fig. \T5\ are plotted against the stretching parameter. 
There is a great deal of scatter in both figures. There is however an observable 
trend, the fragilities appearing to decrease as the stretching increases. The relative 
effect of density (over that of temperature) on the slowing down of the relaxation is 
characterized by the term <^pTg^^^^^ = mp/ntp — 1. In figure [T6l we show the ratio 
mp /ntp as a function of I3kww- Clearly, no correlation is found between this ratio and 
the stretching. 
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Polyvinylchloride 


160, 191 


mill 


140 




0.25 


||62|| 


Polyvinylacetate 


130, 95, 78 


Eiisziiia 


130, 61, 52 


iiisTiiia 


0,43 


|62| 


Polystyrene 


77, 139 




55 


|64| 


0.35 


|62| 


Polymethylacrylate 


102, 122, 102 


BMIisUll 


80, 94 


iiiii 


0.41 


||62|| 



Table 1. Fragilites and KWW stretching exponents of molecular liquids and polymers. The * indicates that the value is not given in the 
corresponding reference but is calculated from the data therein. The following abbreviations are used for the names of the liquids, PC = 
Propylene Carbonate, BMPC = l,l'-bis(p-methoxyphenyl)cyclohexane, BMMPC = l,l'-di(4-methoxy-5-methylphenyl)cyclohexane, KDE 
= cresolphtalein-dimethyl-ether, DEP A = diglycidylether of bisphenol A, and DHIQ = Decahydroisoquinoline. 
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The correlation between stretching and fragility is not strikingly different 
in figures [14] and [151 However, both on theoretical ground (focusing on the 
intrinsic effect of temperature) and on phenomenological one (isochoric fragility and 
stretching do not appear to vary as one changes pressure along an isochrone), our 
contention is that one should prefer using the isochoric fragility. 

In the above we have considered only fragility and stretching at the conventional 
glass transition temperature, that is around = 100 s. However, we have pointed 
out in the Introduction that both the steepness index characterizing fragility and the 
stretching parameter depend on the relaxation time. Although still debated, there 
seems to be a qualitative trend toward a decrease of the stretching (an increase 
in (Sfcww) arid of the steepness index as the relaxation time decreases and one 
approaches the "normal" liquid regime. It would certainly be very valuable to obtain 
more data in order to study how the correlation between fragility and stretching 
evolves as a function of the relaxation time. 
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Figure 14. Isobaric fragility as a function of stretching parameter. Diamonds: 
molecular liquids, circles: polymers. See table l6.2| for numerical values and references. 
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Figure 15. Isochoric fragility nip as a function of stretching parameter. Diamonds: 
molecular liquids, circles: polymers. See table l6.2| for numerical values and references. 
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Figure 16. Ratio between isochoric and isobaric fragility as a function of stretching 
parameter. Diamonds: molecular liquids, circles: polymers. See table 16.21 for 
numerical values and references. 



7. Conclusion 

In this article we have stressed the constraints that one should put on the search 
for (meaningful) empirical correlations between the fragility of a glassformer, which 
characterizes the temperature dependence of the slowing down, and other dynamic 
or thermodynamic properties. Among such constraints is the check that the proposed 
correlations, often established at Tg and at atmospheric pressure, are robust when 
one changes the reference relaxation time (in place of the characteristic of Tg) as well 
as when one varies the pressure under isochronic conditions. Important also is the 
fact that fragility depends on the thermodynamic path considered (constant pressure 
versus constant density) and that, contrary to the isobaric fragility, the isochoric one 
appears as an intrinsic property of the glassformer, characterizing the pure effect of 
temperature. 

We have reported dielectric relaxation spectra under pressure for two molecular 
liquids, m-toluidine and DBF. We have combined these data with the available 
thermodynamic data and analyzed the respective effect of density and temperature 
on the dynamics. Our results are consistent with a general picture in which the 
isochoric fragility is constant on an isochrone. The shape of the relaxation function, 
as e.g. expressed by the stretching parameter I5kwW' has also been found constant 
along isochrones. 

We have finally discussed the possible correlation between fragility and 
stretching, suggesting that a meaningful correlation is to be looked for between 
stretching and isochoric fragility, as both seem to be constant under isochronic 
conditions and thereby reflect the intrinsic effect of temperature. On the practical 
side, the correlation is however no stronger with the isochoric fragility than with 
the isobaric one. One top of large error bars that may be present and that we have 
addressed in some detail, this reflects the fact that correlations are rather statistical in 
nature, emerging from a comparison of a large number of glassformers, rather than 
one-to-one correspondences between properties of the materials. . 
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Appendix A. Details on the density calculation 

The pressure and temperature dependences of the density are of course a crucial 
input to the scaling shown in section I4.1[ In order to evaluate the effect of the 
extrapolations we have performed, we focus on the scaling for the high-pressure 
room-temperature data of Paluch and the data at atmospheric pressure, because the 
extrapolation of the density is smallest in these cases. The discrepancies seen in figure 
in could be accounted for, if the density at high pressure and room temperature were 
higher than what we have estimated or if the density at low-temperature were lower 
than what we have estimated. The high-density dynamical data are taken at room 
temperature. The experimental density data are also taken at room temperature and 
they are only extrapolated above 1.2 GPa. If the actual density is higher than what we 
have estimated, then it means that the compressibility is larger than what we taken. 
However, the compressibility at 1.2 GPa is already in the high-pressure domain 
where it is very low and almost pressure independent (it is slightly decreasing with 
increasing pressure). The most conservative estimate we could make is to keep the 
compressibility constant for pressures above the last experimental point at 1.2 GPa. 
Such an approach changes the ratio p^-^/T by less than one percent, and, therefore, 
can not account for the discrepancy seen in figure |5l An alternative explanation 
would be that the actual low temperature density is higher than we have estimated, 
meaning that we have overestimated the expansion coefficient ap. This latter has 
been calculated at two different high temperatures based on the data in reference ||28|. 
This leads to a slight decrease in expansion coefficient with decreasing temperature. 
If the expansion coefficient is to be smaller than the estimate from this temperature 
dependence, then it would mean that the temperature dependence of the expansion 
coefficient should increase as temperature decreases. This is the opposite of what 
is seen in real liquids, where ocp at atmospheric pressure tends to a constant at low 
temperatures Il30ll . It is actually most common to assume that the ap of molecular 
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liquids is constant below room temperature (e.g. ref. IITSi ). This type of assumption 
would enhance the discrepancy in figure |5l We therefore conclude that the absence 
of collapse of the high-pressure data in Fig. fusing a simple power law form for e{p) 
cannot be explained by errors made in the estimating the PVT data. 

Appendix B. Characterizing the spectral shape 

In the following, we shortly review the procedures commonly used to characterize 
the shape of the relaxation spectrum of viscous liquids and test different descriptions 
on one of our spectra. We more specifically look at schemes for converting one type 
of description to another. This analysis is important for the present work because we 
compile literature data in section [6^ in order to look at possible general connections 
between relaxation shape and temperature dependence of the relaxation time. 

The (normalized) Kohlrausch- William- Watts function or stretched exponential, 
<pKWw{i) = sxp — (^L^^'^^^ ^ leads to a loss peak in the frequency domain that is 
given by the one-side Fourier transform 

cPffKWwicv) = /f -^^^^ sin{cut)dt (B.l) 

The low-frequency behavior of this function is always a power law with 
exponent 1. The high frequency behavior is a power law with exponent —I3kww 
||65|. I3kww is the only parameter describing the shape of the relaxation function. 
Hence it controls both the exponent of the high frequency power law and the width 
of the relaxation function. 

The Havriliak-Negami (HN) function, 

1 

(Phn{co) = TT—J-. 7-T- , (B.2) 

[1 + [icvrHNrr 

gives a power law with exponent (— cty) in the high-frequency limit and a power law 
of exponent ex in the low frequency-limit of its imaginary part. 

The HN function reduces to Cole-Davidson (CD) one when <x = 1. (In the 
case of a CD function we follow the convention and refer to the y above as ^cd) 
The CD spectrum has the same general characteristics as the KWW one: a high- 
frequency power law with exponent given by ^qd ^rid a low-frequency power law 
with exponent one. However, the shape of the two functions is not the same. The CD 
function is narrower for a given high frequency exponent (given ^) than the KWW 
function. The best overall correspondence between the CD-function and the KWW 
function has been determined by Lindsey and Patterson ||65| . 

No good correspondence exists in general between the HN and the KWW 
functions. First of all because the former involves two adjustable shape parameters 
and the latter only one (plus in both cases a parameter for the intensity and one for 
the time scale). The KWW function always has a slope of one at low frequencies 
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while the HN function has a generally nontrivial a. Alvarez et al. ||66ll numerically 
found that the two functions can nonetheless be put in correspondence by fixing the 
relation between the two HN parameters 7=1— 0.812(1 — a)^-^^'^ and choosing 
I^KWW = ((xy)^^/^-'^^\ This restricted version of the HN function is sometimes 



referred to as the A AC function |[67|. The shape is described by one parameter. 
However, it is clear that this function cannot correspond to the KWW function in 
the frequency range where the loss can be described by power laws, as it was also 
noted by Gomez and Alegria ||67| . The A AC function inherits the behavior of the HN 
function; as a result it has a nontrivial exponent a at low frequencies and an exponent 
—ocy at high frequencies, while the associated KWW function has exponents one and 
—I^KWW = — (ay)*-^/^'^^^ at low and high frequencies, respectively. 

Another approach is to describe the dielectric spectrum by a distribution of 
Debye relaxations 

/•oo I 

{e{cv)-eoo)/Ae= / D(/nT)— dlnr, (B.3) 

and to fit the shape of the distribution D{lnr) rather than the spectral shape directly. 
The following form has been suggested for the distribution function ||68B , 

D{lnr) = Nexp{-l3/a{r/ror){r/ro)ni+(—Y ), (B.4) 



where N is a normalization factor. The function above is known as the extended 
generalized gamma distribution, GGE. The last term (and the parameters y and cr) 
describes a high-frequency wing, corresponding to a change from one power law 
behavior {-(5) to another (-y). This term can therefore be omitted if no wing is 
observed in the spectrum. This results in a simpler distribution; the generalized 
gamma distribution (GG) whose shape is described by two parameters: (x determines 
the width and /3 gives the exponent of the high-frequency power law. The low 
frequency is always a power law with exponent one. 

Finally, it is possible to describe the spectra phenomenologically in terms of the 
full width at half maximum, usually normalized to the full width at half maximum 
of a Debye peak |j23 (W/Wd, with Wjj = 1.14 decade), and by the exponent of 
the power law describing the high-frequency side. The power law exponent is not 
always well defined, as there can be a high-frequency wing or a secondary process 
appearing at high-frequencies. Olsen ||69l et al. therefore suggest to characterize the 
alpha peak by the minimal slope found in a double logarithmic plot of the dielectric 
loss as a function of frequency. Note that this phenomenological description requires 
two parameters to describe the shape, while the commonly used CD and the KWW 
functions use only one parameter to describe the spectrum. 

In figure |BT] we show one of the dielectric spectra of m-Toluidtne along with fits 
to the functions described above. The minimal slope is —0.44 and W /Wjj = 1.56. 
The best fits to the different functions are displayed in figure IBli The CD-fit gives 
l3cD = 0.42, which with the Ltndsey-Patterson scheme [65l corresponds to (Bfcww ~ 
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0.55. The direct fit with the Fourier transform of the KWW gives iS^ww = 0.57. The 
best AAC fit gives a = 0.85 leading to y = 0.61 and I3kww ~ (ra;)V^-23 = 0.59. 
This shows that both the Patterson and the AAC approximations reasonably well 
reproduce the ^kww value found from using KWW directly. Another point worth 
noticing is that the I3kww value does not correspond to the actual high-frequency 
slope. This is because the overall agreement between the fit and the data is much 
more governed by the width of the relaxation function than by its high-frequency 
slope, as it is also clearly seen for the KWW fit in figure Ell Note that the AAC 
approximation for the relation between the HN parameters and I3kww only holds 
when the HN parameters are fixed according to y = 1 — 0.812(1 — a)^'^^'^ . The 
original HN function has two adjustable parameters to describe the shape. The best 
HN fit gives a = 0.95, and y = 0.46. The Gamma distribution which also has two 
free parameters gives <% = 40 and (5 = 0.49. Finally we have fitted with the GGE 
using the constraint ^ = 3y (see reference llTOl ), meaning that the function has 3 
free parameters to discribe the shape, the values being a = 40, (5 = 0.7 a = 53 
and y = /3/3 = 0.23. It is not surprising that the GGE with 3 free parameters gives 
by far the best fit. However it is also striking that the CD with only one parameter 
describing the shape gives a good fit over the whole peak, whereas this is not true for 
the KWW nor for the AAC. 

From the above we conclude the CD-function gives a good description of the 
shape of the relaxation using only one parameter to describe the shape. We therefore 
use this function to fit our data. The KWW exponent, ISkv/w^ does not give a proper 
measure of the high frequency slope, but that it does give a reasonable one-parameter 
measure of the overall shape of the dispersion. The KWW function is moreover the 
function most commonly used in literature, which is the main reason for using it in 
the discussion (section (6]). 
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Figure Bl. Log-log plot of the dielectric loss of m-toluidine at T=216.4K and 122MPa 
along with best fits to several common functional forms. Figure a) show the fitting 
functions from below and up; CD, KWW, A AC. Figure a) shows from below and 
up; HN, Gamma distribution. Generalized gamma distribution. CD, KWW and A AC 
have 1 parameter characterizing the shape, HN and Gamma have 2, and Generalized 
gamma has been fitted using 3 adjustable parameters. The dashed line shows the 
Gamma distribution corresponding to the generalized gamma distribution. The 
curves are displaced along the y axis by regular amounts. 



